Neutron-proton bremsstrahlung from intermediate energy heavy-ion reactions as a 

probe of the nuclear symmetry energy? 
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Hard photons from neutron-proton bremsstrahlung in intermediate energy heavy-ion reactions are 
examined as a potential probe of the nuclear symmetry energy within a transport model. Effects of 
the symmetry energy on the yields and spectra of hard photons are found to be generally smaller 
than those due to the currently existing uncertainties of both the in-medium nucleon-nucleon cross 
sections and the photon production probability in the elementary process pn — > pw^. Very inter- 
estingly, nevertheless, the ratio of hard photon spectra -Ri/2(7) from two reactions using isotopes 
of the same element is not only approximately independent of these uncertainties but also quite 
sensitive to the symmetry energy. For the head-on reactions of ^'^^Sn -i-i^"* Sn and Sn -j-^^^ Sn at 
Ebeam/A = 50 McV, for example, the R\/2{l) displays a rise up to 15% when the symmetry energy 
is reduced by about 20% at p = 1.3po which is the maximum density reached in these reactions. 

PACS numbers: 25.70.-z, 24.10.Lx, 13.85.Qk 
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The density dependence of the nuclear symme- 
try energy Esym{p) is important for understanding 
many interesting questions in both nuclear physics and 
astrophysics [H In particular, the reaction dynamics 
and many observables of heavy-ion reactions are strongly 
influenced directly by the corresponding nuclear symme- 
try potential d, U, H]. In turn, these observables can be 
used as effective probes of the symmetry energy Egym {p) ■ 
Indeed, considerable progress has been made recently in 
constraining the Esymi^) at sub- normal densities using 
heavy-ion reactions^, 0, H, Q. However, the behavior 
of the nuclear symmetry energy at supranormal densities 
remains rather illusive. A quick review of the current sit- 
uation of the field reveals that essentially all known 
probes of the Egym (p) in heavy-ion reactions are hadronic 
in nature. These probes, especially if used for study- 
ing the high density behavior of the Esymip), inevitably 
suffer from distortions due to the strong interactions in 
the final state. Although selections of some especially 
delicate observables in certain kinematic/geometrical re- 
gions, such as, the neutron/proton ratio of squeezed- 
out nucleons perpendicular to the reaction plane are 
promising in reducing effects of the final state interaction, 
ideally one would like to have more clean ways to study 
the symmetry energy especially at supranormal densi- 
ties. In this regard, it is interesting to note that the 
parity-violating electron scattering has been proposed 
to measure more precisely the size of the neutron-skin 
in ^°*Pb[l3|. The latter measured previously in experi- 
ments using hadronic probes was shown recently in sev- 
eral calculations to be proportional to the slope of the 
symmetry ene rgy at and/or slightly below the normal 
density [I [H, M, El, El, E3, El, El- Since only elec- 
tromagnetic interactions are involved in electron scatter- 
ings the proposed experiment PREX is theoretically the 
most clean probe known so far for probing the low den- 



sity behavior of the symmetry energy [20|, |2l|]. However, 
practically it is very challenging since the probability of 
obtaining parity-violating electron scattering events in 
the experiment is extremely small. The long waited data 
is hopefully coming soon. Similar to electrons, photons 
interact with nucleons only electromagnetically. Once 
produced they escape almost freely from the nuclear en- 
vironment in nuclear reactions. We notice that soft pho- 
tons from giant dipole resonances in heavy-ion reactions 
have been shown within a semiclassical molecular dynam- 
ics model to be quite sensitive to the symmetry potential 
term in the nucleon-nucleon intcraction[2^. Can one use 
hard photons from intermediate energy heavy-ion reac- 
tions to extract information about the Esym (p) especially 
at supranormal densities? This is a question that one of 
the present authors was repeatedly asked recently by sev- 
eral experimentalists on several occasionsjl^. To answer 
this important question, we report here results of the first 
exploratory study on using hard photons from neutron- 
proton bremsstrahlung in intermediate energy heavy-ion 
reactions as a probe of the Esym{p)- 

Hard photon production in heavy-ion reactions at 
beam energies between about 10 and 200 McV/A had 
been extensively studied both experimentally and the- 
oretically during the last two decades, see, e.g., refs.0, 
[H,!!!] for a comprehensive review. Indeed, very interest- 
ing physics has been obtained from analyzing data taken 
by several experimental collaborations. For instance, the 
TAPS collaboration carried out a series of comprehensive 
measurements at various experimental facilities (GSI, 
GANIL, KVI) studying in detail the properties (energy 
spectra, angular distributions, total photon multiplicities 
di-photon correlation functions, etc), of hard photons in 
a large variety of nucleus- nucleus systems in the range 
of energies spanning Eiab « 20 — 200 MeV/nucleon|27l| . 
They used those bremsstrahlung photons as a tool to 
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study the nuclear caloric curve, the dynamics of nuclcon- 
nucleon interactions, as well as the time-evolution of 
the reaction process before nuclear break-up [1^. The- 
oretically, it was concluded that the neutron-proton 
bremsstrahlungs in the early stage of the reaction are 
the main source of high energy 7 rays. Within the cas- 
cade and Boltzmann-Uehling-Uhlenbeck (BUU) trans- 
port models it was demonstrated clearly that the hard 
photons can be used to probe the reaction dynamics lead- 
ing to the formation of dense matter [2^ [30, HH, [H, [s^l- 
However, effects of the nuclear Equation of State (EOS) 
on the hard photon production was found small [sij. 
While these reaction models were able to reproduce all 
qualitative features of the experimental data, normally 
the quantitative agreement is within about a factor of 
2. One of the major uncertainties is the input elemen- 



tary pn 



arobabilit'v 



which is still rather model 



dependent |35l. l36l. l37l. l38l. l39j . It was noticed earlier that 
the few existing data for the pn — > pn"/ process can be 
described reasonably well by the available models usually 
within a factor of 2[2y| ■ Looking forward enthusiastically, 
we mention here that the very recent systematic measure- 
ments of the pn — > pn"f cross sections with neutron beams 
up to 700 MeV at Los Alamos have the potential to im- 
prove the situation significantly in the near future [40l|. 

Since the photon production probability is so small, 
i.e., only one in roughly a thousand nucleon-nucleon 
collisions produces a photon, a perturbative approach 
has been used in all dynamical calculations of pho- 
ton production in heavy-ion reactions at intermediate 
energies p3.[26j. In this approach, one calculates the pho- 
ton production as a probability at each proton-neutron 
collision and then sum over all such collisions over the 
entire history of the reaction. As discussed in detail 
earlier in ref.jl^, the cross section for neutron-proton 
bremsstrahlung in the long-wavelength limit separates 
into a product of the elastic np scattering cross section 
and a 7-production probability. The probability is of- 
ten taken from the semiclassical hard sphere collision 
modeljli, m, ill. The double differential probabiUty, 
ignoring the Pauli exclusion in the final state, is given by 



We notice that other expressions derived theoretically in- 
volving more quantum-mechanical effects exist in the lit- 
erature, see, e.g., [H, [H, [13, [H, [1^ . Without passing 
any judgement on these theories, to simply evaluate influ- 
ences of the elementary pn — s- pn^/ probability on photon 
production in heavy-ion reactions, for a comparison we 
thus also use the prediction of the one boson exchange 
model by Gan et al.[38j 
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2.1 X 10" 
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where y = Sj/Emax, a = 0.7319 - 0.5898/3^, and E^nax is 
the energy available in the center of mass of the colliding 
proton-neutron pairs. 
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FIG. 1: (Color online) The single differential probability as a 
function of proton kinetic energy in the proton-neutron center 
of mass frame for the production of photons at energies of 
50, 150 and 350 MeV, respectively. The lines with higher 
values are results calculated with the semi-classical Eq. ^ 
while the ones with lower values are obtained by using the 
quantum-mechanical Eq. ^ . 
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^xl(3sin^Mf + 2/5|) 

6.16 X 10"^ X —{3sm^ 9^(3f + 2(3j), (1) 



where 6j is the angle between the incident proton direc- 
tion and the emission direction of photon; and Pi and 
(3f are the initial and final velocities of the proton in 
the proton-neutron center of mass frame. The above 
equation was obtained from modifying the original semi- 
classical Jackson formula 4l| to allow for energy conser- 
vation in the 7-production process [Sll [3^ . Integrating 
Eq. ^ over the photon emission angle, one obtains the 
single differential probability 
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The single differential probability and from the 
two models arc shown in Fig. [T] as a function of pro- 
ton kinetic energy in the proton-neutron center of mass 
frame for the production of photons at energies of 50, 150 
and 350 MeV, respectively. It is seen that the two mod- 
els give quite similar but quantitatively different results 
especially near the kinematic limit where the is signif- 
icantly higher than the p!^, as noticed already in ref . [sst . 
Moreover, the very small magnitude of the photon pro- 
duction probability shown here justifies the use of the 
perturbation method. 

In most of the previous calculations for hard photon 
production using transport models, a constant nucleon- 
nucleon cross section of 30mb to 40mb was normally 
used [13, HH, [13, [13- In this study, we use the 
IBUU04 transport model [l^l where an isospin-dependent 
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in-mcdium nuclcon-nuclcon (NN) cross section 



^medium ^/ree / A*JVJV \2 

fJ-NN 



(4) 



was implemented [8|]. In the above, the /i^jv ^^'^ I^nn 
are the in-medium and free-space reduced NN mass. Be- 
cause of the neutron-proton effective mass sphtting due 
to the momentum dependence of the isovector potential, 
the scaline; factor (^i««-)2 induces ^ significant modifica- 
tion to the relative cross sections of np, nn and pp as 
discussed in detail in ref.0]. The energy and isospin de- 
pendent free-space NN cross section ajj^ are taken from 
the experimental data. 
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FIG. 2: (Color online) Density dependent nuclear symmetry 
energy 

Another important input to the transport model is the 
mean field. In the IBUU04 we use the momentum- and 
isospin-dcpcndcnt single nuclcon potential (MDI) given 
in ref . [431] . In this interaction a parameter x was intro- 
duced to vary the density dependence of the nuclear sym- 
metry energy while keeping other properties of the EOS 
fixed. Fig. [5] shows the density dependent symmetry en- 
ergy with X — 1, X = and x — —I. Available ex- 
perimental data on isospin diffusion Q and isoscalingQ 
have allowed us to constrain the symmetry energy to be 
between the curves with a; = and x = — 1 at subsatura- 
tion densities 0, Q. At high densities, however, there is 
so far no experimental constraint available. One of our 
major motivations here is to examine whether hard pho- 
tons can be used to constrain the symmetry energy at 
supranormal densities as we were asked by the interested 
experimentalists j23j . 

While we are not aiming at reproducing any data in 
this exploratory work, it is necessary to first gauge the 
model by comparing with the available data. Shown in 
Fig. [3] are the calculations with both p'^ and p'^ and the 
experimental data for the inclusive cross section of hard 
photon production in the reaction of ^^C-|-^^C[2^ [il ]. 
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FIG. 3: (color online) Beam energy dependence of the inclu- 
sive photon production cross sections in ^^C-|-^^C collisions. 
The solid symbols stand for experimental data[2^. |43|. (The 
squares are for 50 MeV < £7 < 100 MeV, circles for 100 MeV 
< e.^ < 150 MeV and triangles for > 150 MeV). The solid 
lines are calculated using the and the dashed ones using 
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The calculations are done with a; = 0. It is seen that 
both calculations are in reasonable agreement qualita- 
tively with the experimental data except for the very 
energetic photons. Quantitatively, the agreement is at 
about the same level as p revious calculations by others in 
the literatureflO,!!!, [sl]. We notice that the uncertainty 
in the elementary pn — > pwy probability leads to an ap- 
preciable effect on the inclusive 7-production in heavy- 
ion reactions. As one expects, the larger value of p° 
from the semi-classical picture gives significantly higher 
7-production cross section in heavy- ion reactions. In fact, 
this effect is larger than the effects of the symmetry en- 
ergy obtained by varying the x parameter from a; = to 
X = —1 or a; = 1. It is thus a really very challenging task 
to extract useful information about the symmetry energy 
from photon production in heavy-ion reactions given the 
currently existing uncertainties associated with the ele- 
mentary probability. Nevertheless, we are very hopeful 
that there are ways to overcome this difficulty. As an 
example, we will actually discuss later in this paper one 
possible way to do this. 

We now turn to discussing effects of the symmetry en- 
ergy on photon production in more detail. As an exam- 
ple, we consider head-on collisions of ^''^Sn-l-^^^Sn at an 
incident energy of 50 MeV/ A using the symmetry ener- 
gies of a; = 1 and —1. For this discussion, it is sufficient 
to use only the p^. First, we examine the evolution of 
the central density in Fig. HI The highest central density 
reached is about l.Spo. At this density the symmetry 
energy with a: = — 1 is higher by about 20% as shown 
in Fig. [21 The stiffcr symmetry energy with a; = — 1 
leads to a slightly lower central density. The supranor- 
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FIG. 4: (Color online) The central density reached in the 
head-on reactions ^''^Sn+^'^^Sn at 50 MeV/A with the sym- 
metry energies of x=l, -1 using the p^. The horizontal line 
stands for the normal nuclear density. 
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FIG. 5: (Color online) Time evolution of the multiplicity of 
hard photons with 50 MeV < < 125 MeV in the head- 
on collisions of ^"^■^Sn-l-^'^'^Sn at a beam energy of 50 MeV/A 
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mal density phase lasts from about 10 to 35 fm/c. It is 
in this phase, as shown in Fig. [5l that most of the hard 
photons with 50 MeV < e-y < 125 MeV are produced. 
This observation is consistent with conclusions of the ear- 
lier studies [1^, [13, [13, [13- It is seen that the softer 
symmetry energy with x = 1 produce just a little more 
hard photons. To investigate effects of the in-medium 
NN cross sections, we also carried out calculations using 
the free-space NN cross sections but with the same 
for the case of x = 1. It is seen that the free-space cross 
sections leads to significant higher hard photon produc- 
tion. This is simply because the in-medium NN cross 
sections are smaller than the free-space onces[l]. More- 
over, we notice that the change of the NN cross sections 



leads to a much larger effect than the symmetry energy. 
But overall, both effects are very small compared to that 
due to the uncertainty of the elementary 7-production 
probability. Nevertheless, to put the comparison into the 
proper context we stress that the symmetry energy and 
the NN cross sections arc only varied, respectively, by at 
most 20% and 50% in the reaction considered by chang- 
ing the X parameter between 1 and —1 and the cross sec- 
tions between the free-space and in-medium ones. While 
the current uncertainty of the 7-production probability 
is significantly larger. 
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FIG. 6: (Color online) The spectra ratio of hard photons in 
the reactions of ^^"^ Sn +^^'* Sn and ^^^Sn -|-^^^ Sn reactions 
at a beam energy of 50 MeV/ A with the symmetry energies 
of x=l, x=-l. 

Given all of the uncertainties discussed above and the 
fact that the cross section for hard photon production 
in intermediate energy heavy-ion collisions is very small, 
is there anyway one can reduce these uncertainties nat- 
urally and thus see more cleanly the effects of the sym- 
metry energy? Like in many experiments searching for 
minute but interesting effects, ratios of two reactions can 
often reduce not only the systematic errors but also some 
"unwanted" effects. At least theoretically, within the per- 
turbative approach adopted here, the uncertainty of the 
7-production probability should get almost completely 
cancelled out in the ratio of photons from two reactions. 
We thus propose to measure experimentally the spectra 
ratio Ri/2{l) of hard photons from the head-on reactions 
of i32Sn+i24Sn and "^Sn+ii^Sn, i.e., 



^1/2(7) 



^(132 



124 



Sn) 



^(ii25„+ii2 Sn) 



(5) 



Depending on the relative number of neutron-proton 
scatterings in the two reactions, uncertainties due to 
the NN cross sections can also get significantly reduced. 
Shown in Fig. [5] is the Ri/2{l) calculated for four cases 
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using both the p° and First of all, it is seen clearly 
that the full calculations with the p° and and the 
in-medium NN cross sections indeed lead to about the 
same i?i/2(7) within statistical errors as expected. It is 
also clearly seen that effect of the in-medium NN cross 
sections get almost completely cancelled out. These ob- 
servations thus verify numerically the advantage of using 
the Ri/2{l) as a robust probe of the symmetry energy es- 
sentially free of the uncertainties associated with both the 
elementary photon production and the NN cross sections. 
Most interestingly, comparing the calculations with a; = 1 
and X ~ ~1 both using the p^, it is clearly seen that the 
^1/2(7) remains sensitive to the symmetry energy espe- 
cially for very energetic photons. Again, the symmetry 
energy is varied by at most 20% in the reaction consid- 
ered by varying the parameter x from 1 ro —1. Thus, 
the approximately 15% maximum change in the spec- 
tra ratio represents a relatively significant sensitivity. To 
put things in perspective, we notice that this sensitivity 
to the symmetry energy is actually at about the same 
level as most hadronic probes including the K'^/K^ ra- 
tio. The latter using strange particles is considered as 
among the most clean hadronic probes of the symmetry 
energy [il]. It shows about a 15% change while the sym- 
metry energy is changed by at least 50% at the density 
reached in heavy-ion reactions near the kaon production 
threshold [45}. Obviously, compared to the / K'^ ra- 
tio the hard photon production is an even more sensitive 
and clean observable. Of course, we also notice that while 
photons are completely free from final state strong inter- 
actions, besides cosmic-radiation background, one needs 
to consider photons from tt" and fragment decays in the 
data analysis [ij. 

Naturally, one may think about reactions at higher 
beam energies to reach higher densities and thus to ex- 
plore the behaviors of the symmetry energy there. How- 
ever, other channels for hard photon productions may 
start playing more important roles at higher beam en- 
ergies. Moreover, the reaction dynamics will be then 
dominated by nucleon-nucleon collisions rather than the 
nuclear mean-field. Effects of the symmetry energy on 
photons are then expected to become smaller. This is 
mainly because compared to nucleons that are directly 
influenced by the symmetry potential, the hard photons 
are affected by the symmetry potential only indirectly 



through the momentum distributions and the densities 
of the colliding proton-neutron pairs. This is also why 
the hard photons were found not so sensitive to the nu- 
clear equation of state in an early study [13]. Only at 
intermediate energies both the mean-filed and the NN 
collisions play about equally important roles in the reac- 
tion dynamics. In fact, we also calculated the Ri/2{l) for 
the same reactions but at a beam energy of 400 MeV/A. 
We find a much small symmetry energy effect. 

In summary, we have carried out an exploratory study 
about effects of the symmetry energy on the production 
of hard photons from intermediate energy heavy-ion reac- 
tions using a perturbative approach within the IBUU04 
transport model. Effects of the symmetry energy on the 
yields and spectra of hard photons from individual reac- 
tions are generally smaller than those due to the exist- 
ing uncertainties of both the in-medium nucleon-nucleon 
cross sections and the elementary pn — > pw^ probabil- 
ity. Very interestingly, however, the ratio of hard photon 
spectra Ri/2{l) from two reactions using isotopes of the 
same element is not only approximately independent of 
these uncertainties but also quite sensitive to the symme- 
try energy. The sensitivity is at about the same level as 
most hadronic probes. Compared to the K^/K~^ ratio in 
heavy-ion reactions near the kaon production threshold, 
hard photons are completely free of final state strong in- 
teractions and are even more sensitive to the symmetry 
energy. 
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